Complete conditions for entanglement transfer 
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We investigate the conditions to entangle two qubits interacting with local environments driven 
by a continuous-variable correlated field. We find the conditions to transfer the entanglement from 
the driving field to the qubits both in dynamical and steady-state cases. We see how the quantum 
correlations initially present in the driving field play a critical role in the entanglement-transfer 
process. The system we treat is general enough to be adapted to different physical setups. 
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Quantum networks of remote local processors, which 
are interconnected by quantum and classical channels, 
have been investigated to effectively perform quantum 
computation [jj and quantum communication A 
quantum repeater has also been proposed for error- 
tolerant long-haul quantum communication These 
schemes require a reliable channel to entangle remote 
nodes in order to use in later steps of protocols. The 
usage of a light field to implement a quantum channel 
is a natural choice because of its handiness in generating 
and propagating entanglement Q . However, we have wit- 
nessed that all optical network is technically extremely 
challenging Q. On the other hand, static qubits such as 
the hyperfme structure of atoms are easily accessible and 
manipulative by means of external excitations. There- 
fore, it may be an optimal strategy to use an optical quan- 
tum channel to bring quantum correlation to two remote 
sites of static qubits, where the entanglement is subse- 
quently utilized for quantum information processing, ft 
is, thus, evident that the study of entanglement transfer 
from an optical field of a continuous-variable (CV) sys- 
tem to a static qubit system has a primary importance. 

Recently, such entanglement generation on the pair 
of remote qubits has been studied through the indirect 
interaction via projective measurement two- mode 
squeezed driving field 0, Q and a non-Markovian Q and 
a Markovian environments |lfj| . Although these schemes 
successfully demonstrate the situation for entanglement 
generation on remote qubits, the complete physical re- 
quirements for the possible creation of entanglement are 
unknown. In this paper, we investigate the sufficient and 
necessary conditions to induce entanglement on two re- 
mote qubits, by means of their respective linear interac- 
tions with a two-mode driving field. 

From its definition, entanglement between any two sys- 
tems cannot be created by local unitary operations alone. 
Thus, when the driving field is separable, there is no way 
to generate the entanglement between the two qubits. 
The natural least constraint of the entanglement of the 
qubits is the entanglement of the state for the quan- 
tum channel. However it is not clear if entanglement of 
the driving field can always be transferred to the static 
qubits. We study both the dynamical and the steady- 



state cases. 

Master equation - We analyze the dynamics of two 
remote qubit systems sketched in Fig. ^ Each qubit 
of its ground \g) i and excited \e) i states {i = 1, 2) 
interacts with its own local environment of a single- 
mode bosonic system. We will refer to local environ- 
ments by modes a and b. Static qubits such as ions, 
atoms or quantum dot devices are isolated from un- 
controllable real- world environment. A coupling of the 
static qubit with the driving field is thus assumed to be 
through their small local environments which isolate the 
qubits from the uncontrollable environment. We model 
each qubit-bosonic interaction by a resonant Hamiltonian 
(h = 1), H a i + H b2 , where H al = fi al (&f a + a)a^ ) 
(analogously for H^). Here, a and a' are standard 
bosonic operators, a± = (<rf)t = | e) 1 (e? | and f2 a i is 
the Rabi frequency. The dynamics of the two qubits 
is guided by an external broadband two-mode driv- 
ing field (bandwidth Auj ext ). The coupling between 
modes a, b and the external driving field is written as 

V(t) = Ej «j (aQ^ (t)e^* + aW aj (t)e"^* + a «-» &) , 

where Gaj(t) — aje iujt with a,j the annihilation operator 
of the driving field in frequency ujj , is coupled at rate Kj 
to mode a. Here, u> a ,b are the frequencies of modes a and 
b. In the weak coupling limit Kj <C uJ a ,b, Aw elt (Vj), we 
use second-order perturbation theory and the first Born- 
Markov approximation. The evolution of modes a, b 
can, then, be described by a Liouvillian super-operator 
C involving only the second-order moments of quadra- 
ture variables for the driving field in its carrier fre- 
quency uq, which is set to be resonant to the qubit tran- 
sition, i.e., ujq = Lo a = lo . The coupling rate in this 
frequency is denoted by k. The dynamics of the total 
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FIG. 1: Scheme of the system considered. Two individual 
qubits (Qi, Q2), each interacting with their local environ- 
ment, are driven by a quantum correlated field. 
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system can, thus, be described by the master equation 
d t p = -i[H a i + H b2 ,p] + tp = [c o + Cj p, with p the 
density matrix of the qubits+local environments system. 

In order to characterize the master equation, we need 
to consider only the second-order quadrature moments 
matrix for the driving field which can be written by a ma- 
trix M . The matrix elements are M a p = ({x a , xp}) with 
quadrature operator vector x = (qi,pi,q2,p 2 )- Without 
changing the entanglement structure of the driving field, 
a very general real matrix M can be transformed into the 
simple form by local unitary operations |l2| (no matter 
the driving field being Gaussian or non-Gaussian) 

M=f n C ) (1) 

with n = nl2, m = TOI2 (n,m > 0) that describe the 
local properties of each mode and c = diag[c\, C2] that 
accounts for the inter-mode correlations. 

We are interested in the qubit evolution so to elimi- 
nate the modes a, b. This can be done using an adiabatic 
elimination procedure valid in the weak- coupling regime 
k 3> f2 a i,f2f,2- In this case, the dynamics of the modes 
interacting with the driving field is much faster than their 
interaction with the qubits. The qubits see modes a and b 
in a steady state p ss not affected by the qubit-modes dy- 
namics. The adiabatic elimination proceeds by defining 
a projection operator as V p — p ss <8 Tr a },(p) = p ss ® p X2 . 
Here P12 is the density matrix of the qubits. Using the 
property VC Vp = 0, the reduced master equation takes 

the form d t pi2 = Tr ab ^C a J °° e^C (p ss <g> p X2 ) dt\. It 

is straightforward to find that the dynamics of the qubits 
is fully described by the effective Liouvillian 

4 , s 

CePl2= D °P (VaPuOp- - {6 6 a ,px 2 }) , (2) 

a, (3=1 V - / 

with O a — a a (g> 1 for a = 1, 2 and O a — l®a a -2 for a = 
3,4, cti.2 the x and y Pauli operators. The Kossakowski 
matrix is D = {{A, C}, {C\ BU, where A = A 1 *, B = 
and C are 2x2 matrices [lOj . For the driving field 
with its second-order moments as in Eq. we get 

A = 7i ( n i\ B = 72 / m i \ = V7i72 / ci 0\ 
4 \—i nj 4 y-i m) 4 \ c 2 y 

(3) 

We have introduced the effective decay rate 7^2 = 
2f2^ 1 &2 / K , resulting from the adiabatic elimination. The 
map described by Eq. J2J is completely positive (CP) 
iff D > 0. The interaction model we are using does 
not contain phase-damping processes so that in £ e pi2 
the terms depending on the z Pauli operator are absent. 
Otherwise, the Kossakowski matrix we have, describes a 
general Markovian interaction of two qubits with their 
local environments. C is a real matrix because of the 



constraint of local interaction between qubits and their 
respective environments. 

It is possible to characterize the entanglement capa- 
bilities of the environment-mediated interaction of the 
qubits treated here. In what follows we use the entangle- 
ment measure based on negativity of partial transposition 
(NPT), defined by £npt = — 2Ar, where is the neg- 
ative eigenvalue of the partially transposed density ma- 
trix p^ (T2 indicates partial transposition with respect 
to qubit 2) [l^. NPT is a necessary and sufficient con- 
dition for entanglement of a bipartite qubit system |l5j . 

According to ref. a sufficient condition to entan- 
gle the qubits, which follow the dynamics described by 
Eq. ©, is 

(u f Au) (v f B T v) < jutCvj 2 . (4) 

Here u = (cos 20, — i) T and v = (cos2ip, i) T carry infor- 
mation on the generic initial states of qubits, which are 
unitarily rotated by the angles 6 and ip around the z axes 
of their Bloch spheres. We note that while A, B and C 
depend just on the form of the reduced master equation, 
the condition to entangle the qubits depends on their 
initial conditions via the vectors u and v. We use this 
condition now but will assess it later for the steady-state 
entanglement condition. 

So far the treatment has been very general. However, 
to analyze the entanglement condition clearly, we restrict 
ourselves to the case when c% = — c% = c > in Eq. Q 
from now on. In fact, this case covers most of the entan- 
gled CV states, which have been studied, including Gaus- 
sian noisy two-mode squeezed states and beam-splittcd 
two single-mode squeezed states and non- Gaussian en- 
tangled coherent states after local unitary operations. In 
these conditions and for 71 ^ 72, the map in Eq. (J2J) is 
CP iffc 2 < min{(m-l)(n+l),(m + l)(n-l)}. We find 
that, if the qubits are initially in their ground states, the 
sufficient condition 10} for entanglement becomes 

(n- l)(m- I) < c 2 . (5) 

If this condition is satisfied for the quantum channel, 
entanglement is created between two remote qubits for 
some period of time. The uncertainty principle for the 
quantum channel can be written in the following com- 
pact form: M — cr y © <r y > 0. For the case of a Gaussian 
field, if and only if the partial transposition of their den- 
sity matrix violates the uncertainty principle, the field is 
entangled ^3] and this condition reduces to Eq. (JjJJ 14]. 
Therefore, we find that two remote qubits can be entan- 
gled for some periods during its linear interaction with 
local environments if and only if the Gaussian quantum 
channel of \c\\ — \c 2 \ is entangled by appropriately choos- 
ing initial states of the qubits . The second-order mo- 
ment matrix is obtained from a general case by local 
unitary operations of the field. This can be interpreted as 
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transforming the initial qubit states as leaving the quan- 
tum channel in the simple form Q . For the dynamic en- 
tanglement of the qubits, we have to prepare the initial 
states of the qubits carefully. For a non-Gaussian field, 
the uncertainty principle serves only a sufficient condi- 
tion of its entanglement. 

Let us investigate the entanglement condition In 
order to do it, we further assume a special case of n = m 
and 71 = 72 = 7 for a short while as it is not straightfor- 
ward to solve the general dynamic equation J5J. In this 
case, Eq.(jS| becomes n — 1 < c. For the density matrix 
elements pijhk — {ij\p\2\hk) (i,j,h,k — e,g), we find the 
coupled Bloch equations 

Peeee — 7 [ 2il^p eeee ~t~ ^liPegeg ~t~ Pgege) ~t~ CPeegg\ : 
Pegeg = 7 [«1 (1 ~ Pgege) + P 
Pgege = 7 — Pegeg) + P 

Peegg = ~7 [ n 2Peegg ~ c(l/2 

where, n l 



'eeee ^sPegeg CPeegg] 
( ee T^sPgege £Peegg\ 
qege — Pegeg)] , 



k , 



(6) 

k — -^{n — 1) + I and, by hermiticity, p gge e = 
Peegg- All the other matrix elements are decoupled from 
these equations. The normalization condition determines 
Pgggg. Solving Eqs. © , we find the dynamics of the 
entanglement between the two qubits as shown in Fig. El 



(a) 



(b) 



•'npt 



0.8 
0.6 
0.4 
0.2 



& 

0.6 
0.4 
0.2 



Pi 



20 Z= 40 = 60 = 80 



FIG. 2: Dynamics of entanglement for the two qubits. The 
entanglement is plotted as a function of the dimensionless 
interaction time r = yt for n — 2.4 and three different values 
of c: c = 1.58 (solid curve), c = c ss (2.4) = 1.804 (dashed 
curve) and c = \/n 2 — 1 = 2.18 (dot-dashed curve) 
initial state is \gg) 12 (a) while it is |ee} 12 (b) . 
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By inspection of Fig. [3(a), it is apparent that, even 
for n — 1 < c, the long time behavior of the entangle- 
ment function can lead to a separable qubit state. This is 
due to the fact that the condition JSJ for the dynamical 
entanglement-transfer does not give information about 
the steady-state entanglement. The criterion (0} for en- 
tanglement due to the interaction with a Markovian en- 
vironment is, indeed, relative to the creation of quantum 
correlations in an initially separable state. The sufficient 
condition J3J to entangle two qubits comes from a pos- 
itive gradient of £npt at t = 0. This is the case for 
pi2(0) = \gg) 12 (gg\ but not for pi 2 (0) = |ee) 12 (ee|, for 
example. For this initial state, dtSNPT < at t = 0, as 
seen in Fig. [21(b). In fact, if we start with excited states 



for qubits the sufficient condition leads to c > n + 1 , 
which is physically meaningless as the CP condition and 
the uncertainty principle impose c 2 < n 2 — 1 (we remind 
that n > 0). Does it impose that qubits initially in their 
excited states will never be entangled during their dy- 
namics? The answer is 'no'. The condition i§ only 
a sufficient condition and Fig. (b) clearly shows that 
the qubits can be entangled at a later period of the dy- 
namical evolution. A similar result is found for when one 
qubit is prepared in \g) while the other is in |e). 

It is, thus, interesting to investigate the conditions to 
entangle qubits in their steady state. We now lift the 
temporary condition n = m and 71 = 72 while keep- 
ing |ci| = |c 2 | and find the boundary value c ss (n, m) 
of the correlation parameter c beyond which we are 
sure that the qubit steady state is inseparable. To find 
c ss (n,m) we have to look at the asymptotic behavior 
of the entanglement function that can be found solving 
the generalization of Eqs. @ to n ^ m, 71 7^ 72 and 
looking for steady-state solutions. Then, the condition 
limt_>oo £NPr(t, n, m, c)| c=Cas ( n)m ) = fully character- 
izes the boundary value. We find 



(n, to) 



7171 + TO7 2 
2V7172 



/Lt(n,m) — yj v(n, m) 
n 2 m 2 + (to — n) 2 



(7) 



with fi(n, to) = (nm — l) 2 + (nm + 1) — (n — to) 2 and 
v(n,m) — Anm + A(nm — l) 2 — 3(rt — to) 2 . The two 
qubits are entangled at their steady state if and only if 
c > c ss {n,m). 

Cavity quantum electrodynamic (CQED) system - We 
consider a CQED setup to illustrate the conditions for 
efficient entanglement transfer. This model was recently 
suggested by Kraus and Cirac Q to show the possi- 
bility to entangle two identical two-level atoms respec- 
tively trapped in two remote single-mode cavities. The 
cavities are driven by a broadband two-mode squeezed 
state. Here, we show that our general approach gives the 
complete conditions to entangle the atoms. The cavity- 
driving field coupling rates n are taken to be identical 
under the identical cavity assumption. The adiabatic 
elimination described above gives us the reduced atomic 
master equation. The weak-coupling regime is now equiv- 
alent to the bad- cavity limit in which the steady state of 
the cavity is the two-mode squeezed state as in the case 
without atoms in the cavities. Experimentally, k can not 
be taken large at will because the relation Auj ex t 3> K 
represents a constraint for the validity of our treatment. 
Typically, it is Aui ext ~ k/6 ~ 2tt x 12 MHz and val- 
ues as (fi,r)/27r ~ (20, 3.5) MHz allow for the validity 
of the bad-cavity regime and for the squeezed state to 
build up inside the cavities 0. For the CQED sys- 
tem here, the atoms are entangled in the steady state 
if c > c ss = i (^/ (n 2 — l) 2 + n 2 — 1^. This is a severe 
constraint on the properties of the driving field. The ex- 
perimentally available source of squeezed light is, indeed, 
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quite bright, i.e., n large, so that c ss — > \Jn 2 — 1. This re- 
stricts the range of values of c in which the atomic steady 
state is entangled. Note that c 2 = n 2 — 1 is satisfied by a 
pure state of the driving field. 

In the CQED model, the atoms may interact not only 
with the single-mode cavity fields but also with other 
uncontrollable reservoir through atomic spontaneous de- 
cay. Including the atomic spontaneous emission of its 
rate T, the Louivillian remains in the form J2J but with 
71,2, n and c replaced by j[ 2 = (2£l 2 alb2 / + 1/C), 
n' = ml = (nC + 1)/(1 + C) and d' = cC/(l + C) 
(C = 2Q 2 /Tk is the cooperativity With these new 

parameters and for the atoms initially in their ground 
states, the condition to entangle the two dynamic atoms 
still remains as (n — 1) < c for the driving field. We find 
that even with the atomic spontaneous decay, the atoms 
are guaranteed to be entangled for periods of time by 
interaction with the entangled squeezed field jl8j. 

With these new parameters, it is always n' 2 > c' 2 + 1, 
even for a pure drive. Thus, inside the cavities, a pure 
two- mode squeezed environment, where the discussion 
in [j| was centered, is hard to be obtained in a realis- 
tic situation. It is thus worth addressing the effect of the 
purity of the quantum channel in this example. As a mea- 
sure for the purity of the atoms, we take the linearized 
entropy S L = 4/3 (l - Tr 12 {p 2 12 }) that ranges from 
(pure states) to 1 (maximally mixed ones). Only the 
interaction with a pure squeezed environment realizes a 
pure atomic steady state. The dynamics of the linearized 
entropy Sl f° r the atoms initially in their ground states, 
are shown in Fig. |5J The higher is \Jn 2 — c 2 , the more 
mixed the driving field and the higher is (n— c) _1 and the 
more the environments are entangled 14] . It is apparent 
that only a slight departure from the pure state of the 
driving field brings about the atoms extremely mixed as 
shown in Fig. However, as stated before, in a realistic 
situation, V ^ so that it is not possible to get the local 
environments in their pure correlated state. Hence, the 
atomic steady state will always be mixed. When C ^S> 1, 
a nearly pure steady state may be obtained, however, the 
bad cavity limit may not be used and the probability to 
feed the cavities becomes small [lfij |. 

Remarks - In this paper, we investigated the conditions 
to entangle two remote qubits dynamically and in their 
steady state, addressing the case of Markovian interac- 
tion with their local environments, which are driven by a 
quantum channel. We found that the entanglement of the 
Gaussian environments is not only a necessary but also a 
sufficient condition to see entanglement of the qubits for 
some period of their evolution, provided the qubits are 
appropriately prepared at the initial instance. We found 
the boundary value of the correlation parameter of the 
quantum channel only beyond which the steady state of 
the qubits is entangled. 
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FIG. 3: The linearized entropy Sl is plotted against r = 
yt for the pure environments \/n 2 — c 2 = 1 and correlation 
(n — c) _1 = 2.5 (solid line) and for the mixed environments 
\/n 2 — c 2 = 1.4 and (n — c) _1 = 2.5 (dot-dashed line) and 
(n — c) _1 = 3.5 (dashed line). 
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